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Abstract. A 3D simulation of the upper part of the solar convective zone is used 
to derive constraints about the averaged and dynamic properties of solar turbulent 
convection. Theses constraints are then used to compute the acoustic energy supply 
rate P injected into the solar radial oscillations according to the theoretical expression 
in Samadi & Goupil (2001). The result is compared with solar seismic data. 
Assuming, as it is usually done, a gaussian model for the frequency component Xk( v ) 
of the model of turbulence, it is found that the computed P{y) is underestimated com- 
pared with the solar seismic data by a factor ~ 2.5. 

A frequency analysis of the solar simulation shows that the gaussian model indeed does 
not correctly model Xk( u ) i n the frequency range where the acoustic energy injected 
into the solar p-modes is important (y ~ 2 — 4 mHz). One must consider an additional 
non-gaussian component for Xk(v) to reproduce its behavior. Computed values of P 
obtained with this non-gaussian component reproduce better the solar seismic obser- 
vations. This non-gaussian component leads to a Reynolds stress contribution of the 
same order than the one arising from the advection of the turbulent fluctuations of 
entropy by the turbulent motions. 
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1. Introduction 

Turbulent motions in the upper convective zone of solar like stars, generate acoustic energy, 
which in turn is injected at the rate P{y) into p-modes oscillations [y is the frequency of a 
given mode). Solar- type oscillations are therefore meant as stochastically excited by turbu- 
lent convection. This type of excitation process concerns low massive stars with an outer 
convective zone. 

Measurement of the rate of acoustic energy, P{v), injected into solar-like oscillations is 
one of the goals of future space seismic missions COROT (Baglin et al 1998) and Eddington 
(Favata et al 2000). These seismic data will then make it possible to constrain the the- 
ory of excitation and damping of solar-type oscillations, which in turn will provide valuable 
information about the properties of stellar convection zones. 

Theoretical formulations for P{v) offer the advantage of testing separately several proper- 
ties entering in excitation mechanism. Here we consider the formulation by Samadi &; Goupil 
(2001, see section ^] and also Samadi (2001) for a detailed summary). 

Computation of P{v) requires an accurate knowledge of the dynamic properties of turbu- 
lence in stars. Unfortunately, the current observations of the solar granulation cannot provide 
a precise enough description of the turbulent spectrum properties. 

In the present work we use a 3D simulation of the upper part of the solar convective zone 
to determine the dynamic properties of solar turbulence (see section |3]) which are necessary 
to compute P{v)- Our calculations are then compared with the helioseismic constraints from 
the GOLF/SOHO instrument (see section ^]). 



2. A theoretical formulation for P{v) 

The calculation of the rate P(v) at which a given p-mode is excited results from an integration 
over the stellar mass (m) and local integrations over distance (r) and time (t) of the mode 
eigenfunction (£) and the correlation product of the excitation source (< SS >). This 
expression can be written in a schematic form as : 

P(u) oc J dm J drdr £. < SS> .|* (1) 

The excitation source (S) has two identified origines : the turbulent Reynolds stress 
and the advection of the turbulent entropy fluctuations by the turbulent motions. < S S > 
is expressed in terms of the turbulent kinetic energy spectrum E(k, v) and the spectrum of 
the entropy fluctuations E s (k,v) where k is the wavenumber of a given turbulent element. 
Detained expressions for P{v) and < S S > are given in Samadi & Goupil (2001). 

Following Stein (1967), E(k,v) is split into a spatial component E(k) and a frequency 
component Xk(v) as 

E{k,u) = E{k) X k{y) ■ (2) 

The same decomposition is assumed for E s (k,v). 

A gaussian shape for Xk{v) is usually assumed in the calculation of P{v) (e.g. Goldreich 
& Keeley 1977). This is equivalent to suppose that two distant points in the turbulent 
medium are uncorrelated. In section ^] below, the gaussian hypothesis is compared with the 
properties of Xk{ v ) inferred from a 3D simulation of the solar convective zone. Section 4 
discusses consequences of either choice on the excitation rate P{v)- 
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3. Constraints from a 3D simulation of the upper part of the solar convective 
zone. 

We study a 3D simulation of the upper part of the solar convective zone obtained with the 3D 
numerical code developed at the Niels Bohr Institute for Astronomy, Physics and Geophysics 
(Copenhagen, Denmark). Physical assumptions are described in Stein & Nordlund (1998). 

The simulated domain is 3.2 Mm deep and its surface is 6 x 6 Mm 2 . The grid of mesh 
points is 256 x 256 x 163, the total duration 27 mn and the sampling time 30s. 

The simulation data are used to determine the quantities E(k,z), E s (k,z) and Xk( v ,z) 
involved in the theoretical expression for the excitation rate P{y). 

We proceed in two steps (more details will be given in Samadi et al 2002a and Samadi 
et al 2002b): 

• We compute at each layer z the 2D Fourier transform, along horizontal plans, of the 
velocity field u and the entropy s and perform integrations over circles with radius k. 
This provides u(/c, z, t) and s(k, z, t) where k is the wavenumber along the horizontal 
plan. We finally time average each quantities u and s over the time series. This provides 
u(k,z), s(k,z) and then time averaged kinetic energy spectrum E(k) = u 2 (k,z) and 
the time averaged spectrum of the entropy fluctuations E s (k) = s 2 (k,z). 

• At 3 different layers of the simulated domain, we compute the 3D Fourier transform, 
with respect to time and along the horizontal plan, of the velocity field u. We next 
perform integrations over circles with radius k. This yields u(/c, z, v) and therefore 
E(k, v, z) = u 2 (k, v, z) and - using Eq.(2) - Xfc(^, z). 

The dependence of Xk with frequency is plotted at the top of the superadiabatic region 
-where the p- modes excitation is the largest - (Fig. 1), 0.3 Mm deeper (Fig. 2) and 0.6 Mm 
deeper (Fig. 3). 



We compare Xk{ u ) obtained from the simulation analysis with three analytical functions. 
These are the Gaussian function (GF hereafter): 

Xkiy) = -^-=e~^ 2 , (3) 
Vk V 71 " 

the Gaussian plus an Exponential function (GEF hereafter): 

»M(^ _w ^ + £ eH ' M ')- (4) 

and the Gaussian plus a Lorentzian function (GLF hereafter): 

where vy. is the line- width at half maximum of Xk- v k is related to the characteristic 
correlation time-scale of an eddy with wavenumber k as 

Vk = (TTTfe)" 1 (6) 

and Tfc is related to the velocity Uk of an eddy with wavenumber k as 

r k = X(ku k )~ 1 . (7) 
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Figure 1. The solid curve with dots represents Xk{ v ) obtained from the simulation 
at the top of the superadiabatic region (z = 0.04 Mm) and for the wavenumber k 
at which E{k) is maximum. The solid curve represents the Gaussian function (GF, 
Eq. ||), the dashed curve the Gaussian Exponential function (GEF, Eq. ||) and the 
dots-dashed curve the Gaussian Lorentzian function (GLF, Eq. |3|). 
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The velocity Uj- is obtained from the kinetic energy spectrum E(k) (Stein 1967) as 

4= dkE(k). (8) 



The parameter A in Eq. [7| accounts for our lack of precise knowledge of the time correlation 
(or ffc) in stellar conditions. We find however tgat in most part of the excitation region, 
the line width of Xk is satisfactorily reproduced with A = 1. 

At the top of the superadiabatic region, the GF and GLF do not correctly model Xk( u ) 
whereas the GEF is in better agreement (see Fig. 1). However the discrepancies between 
the GF (or the GLF) and the simulation data occur mostly above the solar cut-off frequency 
(y ~ 5.5 mHz). Discrepancies between the GF (or the GLF) and the simulation data have 
then minor consequences for the p-modes excitation in this region. 

This is not the case deeper in the simulation where the largest discrepancies between 
the GF and data occur in the frequency range where the larger amount of acoustic energy is 
injected into the p- modes (y ~ 2 — 4 mHz). The GEF and GLF reproduce better than the GF 
the ^-variation of Xk below the top of the superadiabatic region (z < Mm, see Fig. 2 & 3). 



4. Consequences in terms of p-mode excitation 

We compute P{v) according to Eq. (1): 

• The eigenfunctions (£) and their frequencies (y) are computed with Balmforth's (1992) 
non-adiabatic code for a solar ID mixing-length model based on Gough's (1977) non- 
local time-dependent formulation of convection. 

• The /c-dependency of E(k, z), is modeled as following : 

E(k) oc {k/k ) +1 for k > k > 0.17 k 



E{k) (x (k/k )- 5 / 3 for k > k 



o 



(9) 



where ko = 2ir/(3A, A = aH p is the mixing-length, H p the pressure scale height and 
a the mixing-length parameter. The value of a is this of the ID solar model for con- 
sistency.The value of ko hence of j3 is obtained from the simulation. This analytical 
/c-dependency of E reproduces the global features of E arising from the simulation. 
Same model is considered for E s (k, z). 

E(k,z) and E s (k,z) verifie the normalization conditions: 

JdkE(k,z) = - < u 2 - < u > 2 > (z) 

1 (10) 

dkE s (k,z) = - < s 2 - < s > 2 > (z) 



I 



where < . > denotes time and horizontal average. The total energy contained in E(k, z) 
and E s (k,z) and their depth dependences are then obtained from the simulation ac- 
cording to Eq.(10). 



• For the frequency component Xfc(^)> we assume successively the GF, the GEF and the 
GLF (see section ^]). 
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Figure 4. The curves correspond to computed P{y) in which we assume different 
anaiyticai functions for Xk{ v ) : the GF (solid curve), the GEF (dashed curve) and 
the GLF (dots-dashed curve). The dots represent P(v) derived from the amplitudes 
and line widths of the £ = p- modes measured by GOLF/SOHO instrument and 
kindly provided by F. Baudin (Baudin et al 2003, see also Thiery et al 2000). No 
significant difference for our purpose here are observed between the GOLF/SOHO 
data and Libbrecht (1988) observations. 
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Results are presented in Fig 4. Using the GF leads to a significant over-estimation of 
P(v). This is a consequence of the large discrepancy between the GF and the 3D simulation 
frequency dependence of Xk{ u )- Using the GLF yields a power P{v) which is much larger 
than the observations in particular at high frequency. This is because at the top of the 
superadiabatic region - where the excitation is the largest - the GLF over-estimates Xk(y). 

The best agreement between computed and observed P(v) is found when assuming the 
GEF. This is a direct consequence of the rather good agreement between the GEF and the 
dynamic behavior of the solar turbulence inferred from the simulation (see Sect. 3). 

5. Conclusions 

Our investigation demonstrates the non-gaussian character of the stochastic excitation of solar 
p-modes. Indeed, the gaussian function (GF) as a model for Xkiy) leads to an important over- 
estimation of P{v) whereas the Gaussian Exponential Function (GEF), which decreases more 
slowly with v than the GF does, results in a better agreement between P(y) and the seismic 
solar observations. The maximum value of P{y) is now reproduced without any adjustement 
of free parameters in contrast with previous approaches. 

The non-gaussian character of stochastic excitation causes the Reynolds stress contribu- 
tion to be of the same order as the contribution arising from the advection of the turbulent 
fluctuations of entropy by the turbulent motions (not shown here, see Samadi et al 2002b). 
This last result is in better agreement with recent results by Stein & Norlund (2001) and 
contrasts with previous results by Samadi et al (2001) based on the gaussian assumption for 
Xkiy) an d also with the estimation carried out by Goldreich et al (1994). 
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